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Dimuonium (the bound system of two muons, /x"'"/i~-atom) has not been ob- 
served yet. In this paper we discuss the electromagnetic production of dimuo- 
nium at RHIC and LHC in relativistic heavy ion collisions. The production of 
parastates is analyzed in the equivalent photon approximation. For the treat- 
ment of orthostates, we develop a three photon formalism. We determine the 
production rates at RHIC and LHC with an accuracy of a few percent and 
discuss problems related to the observation of dimuonium. 
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I. INTRODUCTION 



The study of exotic electromagnetic bound systems and their properties is of theoretical 
and experimental interest. The bound /i system (dimuonium, DM) has been subject 
to extensive theoretical investigations As demonstrated in 0, the decay rate of the 

dimuonic system is sensitive to radiative corrections from the so far unexplored time-like 
region of QED. 

Although dimuonium has not been observed yet, a lot of different pathways for its pro- 
duction have been considered. For example, the production of dimuonium in the decay of the 
77-meson [f] —* DM+7) was investigated in Refs. and in Ref. [0] the decay DM+7 
was considered. It has to be mentioned that in decays it is possible to produce only the S = 1 
orthostates of dimuonium. Other calculations were performed for the production of dimuo- 
nium in collisions of charged particles (see Refs. and in collisions of photons with nuclei 

i- 

In this paper we investigate quite a different mechanism, which is based on the availability 
of relativistic heavy ions at high luminosities. Two new large hadron colliders, RHIC and 
LHC, are scheduled to be operative for the next decade. In Table | we list the decisive 
experimental parameters of the new colliders (see Refs. p|-p!0|). We consider here the purely 
electromagnetic production channel 

Ai A2 ^ Ai + DM , 

where the Ai represent relativistic nuclei with nuclear charge numbers Zi, and DM stands 
both for the S = parastates of dimuonium (paradimuonium, PM) and for the S = 1 
orthostates (orthodimuonium, OM). Because the nuclei do not change during the production 
process, they emit the photons coherently. This means that the perturbation parameter 
associated with each photon exchange between the nuclei and and the produced system is 
not a ~ 1/137, but rather Za ~ 0.6 (for Au and Pb). This leads to a very large flux of 
equivalent photons available for the production of exotic particles. 

The C-even PM can be produced in collisions of an even number of virtual photons (two 
photon production mechanism, see the diagram of Fig. |l]a. The C-odd orthostate (OM) can 
only be produced by an odd number of virtual photons, i.e. via bremsstrahlung production 
(one photon. Fig. |l]b) and three photon production (see Fig. |l]c). We consider here mainly 
the production of PM by two photons and the production of OM by three photons. Two 
photon and three photon fusion is the dominating process for the production of parastates 
and orthostates, respectively. The influence of multiphoton processes on the production rate 
is described by the effective perturbation parameter 



fZaAV 



with 



0.04, (1) 



1/A^ = 1/6 {r') , (2) 

where (r^) is the mean square radius of the charge distribution of the nucleus, and the mass 
of the dimuonic atom is m^^ ^ 2 = 211 MeV. Therefore, in all cases under consideration 
the multiphoton processes set limits on the accuracy on the level of 5 %. 
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For our purpose it is sufficient to treat the dimuonia as compound neutral particles. To 
a good approximation their production rate then is proportional to the square of the wave 
function at the origin. In non-relativistic approximation it is only the probability density of 
S states at the origin which does not vanish, 



2 _ a^ml 



\M0)\' = —^. (3) 



The production rate and the lifetime of the dimuonic atoms are both proportional to this 
value. The lifetimes of low lying states are of the order of r ~ 10~^^ s and are summarized 
in Table 0. A brief discussion of the evaluation of the lifetime of parastates is given in 
Appendix 0. Main decay channels are the annihilation processes 

PM ^ 77 , OM ^ e+e~ . (4) 

The rate of atomic transitions from excited S states to lower atomic states is of the same 
order of magnitude {a^ m) as the annihilation decay rate. It results in additional final states 
via atomic decays of excited DM levels which cascade through S — P — S transitions. 
This leads to observable X-ray photons (at least two quanta) having "atomic" energy ~ 
(n'^^ — n^^)/4. Main properties of the various states can be found in Table |I| together 
with 2P paradimuonium which is produced in atomic transitions from 3S and 4S. 

The detection of dimuonium would constitute a continuation of the recent investigations 
of exotic bound systems. Over the past years, experiments on antihydrogen []rT],|T2[, pionium 



||13| and the bound Tifi system |114| , p!5[| have been reported. 

This paper is organized as follows: first we investigate the production of paradimuonium 
in Sec. |I[ We then proceed to orthodimuonium, which is discussed in Sec. |I|. Finally we 



discuss the background in Sec. IV and summarize the results in Sec. 0. 



II. PARADIMUONIUM PRODUCTION 

The production of an S* = parastate of dimuonium by a two photon process is rep- 
resented by the diagram in Fig. |. The diagram is evaluated using the equivalent photon 
approximation in the approach originally presented in Ref. [|T6|. Two nuclei Ai and A2 
with identical charge number Z and atomic mass number A colliding with each other emit 
drii , (z = 1,2) equivalent virtual photons within the energy ranges {uji,uji + dui) and with 
four-momenta denoted as qi. The virtualities of the photons are Qf = —qf. Upon fusion, 
these photons produce a PM-bound state with four-momentum p = qi + q2- Its mass squared 
p2 = _ ^g,^ _|_ g,^^2 jg approximately equal to 4m^. The most important contribution to 
the production process stems from photons with very small virtualities -C m^. To a good 
approximation, the photons move in opposite directions, and we have ~ Aujiuj2- In this 
very region the differential cross section da for the Ai A2 A1A2 + PM process is related 
to the cross section a^^ for the process 77 PM by the equation 

•^CpM = dnidn2 a^^^iW'^) . (5) 



The spectrum of equivalent photons is given by Eq. (D.4) in Ref. |]T6[, which upon omission 
of terms of order uJi/E reads 
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IT UJi \ Qi / Qi Y 



In the calculations below we do not use the exact form factor of the nucleus F^Q"^) but a sim- 
ple approximation. This approximation corresponds to an exponentially decreasing charge 
distribution of the nucleus, whose mean square radius is adjusted to fit the experimental 
value (see Ref. 0, Eq. (B49)). 



i^^r^2^ 1 u a2 0.164 GeV' _ 

According to Eq. (Q), for Pb and Au the parameter A 70 MeV, and for Ca A 118 MeV. 
The approximate form factor enables us to perform some calculations analytically which else 
could only be done numerically. 

It is useful to note that the integral over converges fast at > h?. The decisive 
region of integration is given by the condition Q^^^ < Q"^ ^ (cf. Eq. (|)). Therefore the 
main contribution to the cross section is given by virtual photons with energies 

UJi ~ A7 . (8) 

Because the two photon width of paradimuonium is small in comparison with its mass, 
we can use a 5-approximation for the cross section a^^ (for further details see Eq. (3.24) in 
Ref. [|l^ and Eq. (89.4) in Ref. ||18|). For the l^S'o para ground state, this approximation 



has the form 



After the transformation 



a,,{p') = 2ir'aH{p'-Aml). (9) 



^ ^ 5(/ - 4 = ^ ^ 5{p' - 4 

UJl UO2 UOl p 

1 duJi 
4 oji 

we cast the cross section into the form 

TT^ 5 duJi dni{uji, Ql) ^^2(^2, Ql) 



(10) 



dapM = — a" 



where UJ2 = im?^/uji. Using this formula, we derive the distribution of the produced PM atoms 



2 uoi duoi duj2 

with respect to the energy e and the transverse momentum p^ via the relations 

e = uji^ ^ , p± = qi± + q2±. (11) 

UJl 

It is useful to note that the integral over Q"^ converges fast for > A^. Integrating driiiuj, Qf) 
over Ql, we obtain the equivalent photon spectrum in dependence on the energy, driiiuj). 



dni[uji) -- 



'^a J ujA duJi 

vr HatJ ujr ^ ^ 
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The function 

/(x) = (1 + 2x2) In (l + l) -2 (13) 

drops very quickly at large x in accordance with Eq. (||) (indeed, /(x) < 1/(6 x^) at x > 1). 
Finally, we obtain 

apM = |^G(5), where 6 = ^ (14) 
2m^ A7 



and 



'^min 



-/(x5)/(5/x), (15) 

X 



with X = uji/m^. Because Ui < E and co'iCo'2 = we have x^m = ''tt-^i/E and x^aa; = E/m^. 
However, due to the fast decline of /(x) at x > 1 we can expand these limits up to Xmin = 
and Xmax = 00 in very good approximation. 

Numerical evaluation of the integral in Eq. ([T5|) yields the following result for the total 
production cross sections 

0.15 for RHIC, Au mode , 
(XpM = 10"^° cm^ X { 1.35 for LHC, Pb mode , (16) 

0.0066 for LHC, Ca mode . 

The production cross sections for excited nS states are derived from the above cross section, 
which is obtained for the IS states, with the aid of Eq. (^, 

a{nS) = ^ . (17) 

The summation over n enhances the result of Eq. (|16D by a factor of 

00 

C(3) = E-^ 1-202. (18) 

n=l 

The distribution in energy for the paradimuonium atoms is given by the integrand of Eq. (|15|), 
using the relation Eq. (0). It is shown in Fig. |^. 

We checked that the results depend only weakly on the choice of the form factor. With the 
help of a numerical computer program |19| we took into account the Gaussian form factor 
exp(-QVA2) with A^ = 60MeV (for Pb and Au collisions) and Kg = 100 MeV (for Ca 
collisions) fixed on (r^). We found that this changes the final result presented in Eq. ([T6|) by 
less than one percent. The effect of omitting terms of the order 0{Qf /mjj) in the equivalent 
photon spectrum Eq. (^) is also negligible. The relative contribution of the omitted terms is 
of the order of 
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A2 1 A^^^ 

r]2 = ^ — , and with L = In — = In — — it follows 1/2 ~ (1 -i- 2)% . (19) 

2mj, L 0^ mj, 

The effect of processes involving more than two photons is governed by the parameter p 
from Eq. (0) and we recall that their effect is on the level of roughly 5 %. 

It is instructive to consider additionally the leading logarithmic approximation (LLA) 
for the process. In the LLA, we approximate f{x) by 2 ln(l/x). The restriction min ~ 
corresponds to m^/ (A7) < ui < A7. Therefore 

= -L^ , (20) 
3 



and 



^PM • (21) 

The above result is in good agreement with the old result of ||20| (see also Eq. (2.4) in review 
p!6|). However, for the energies discussed in this paper the LLA does not provide sufficient 
precision. The ratio G^^^/G is 1.5 for Pb at LHC and 2 for Au at RHIC. Hence, the LLA 
gives only a crude estimate for the energies discussed. 

III. ORTHODIMUONIUM PRODUCTION 

Orthodimuonium can be produced by bremsstrahlung (the relevant diagram is depicted 
in Fig. |T]b) and by three photon fusion (see Fig. |I]c). For production processes induced 
by relatively light particles (like e^e'^ or pp collisions) the three photon cross section a^^ 
corresponding to Fig. |l]c is suppressed by a factor compared to the cross section for 
bremsstrahlung production, dbr- By contrast, for heavy ion collisions another parameter 
enters the calculation: the large nuclear mass M. Bremsstrahlung of heavy particles is sup- 
pressed by a factor l/M^, so we obtain 



r2 

5 



dbr oc Z^'a'/M' 

whereas for three photon production there is no such suppression, 

CT3T, OC Z^a^/ml . (22) 

The ratio 

^ 2 f 1/150 for RHIC, Au mode, 

^ < i_ ( = J 1/190 for LHC, Pb mode, 

"^3^ a VM/ y for LHC, Ca mode 

is small. Moreover, a more accurate estimate for CaCa collisions at LHC decreases this ratio 
at least by a factor of three. Therefore, the 3 photon production dominates in relativistic 
heavy ion collisions. In Fig. only one representative diagram for three photon fusion is 
depicted. For a complete analysis, we need to take into account two classes of diagrams, in 
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which the single photon is emitted by either one of the nuclei (see Fig. and Fig. |^b). The 
corresponding cross section, which is proportional to the square of the amplitude for the 
particular processes, is given by 

d<JoM = d(Ta + d(Tb + C^CIinterf = 2 da a (23) 

because the interference term (io"interf disappears after azimuthal averaging. Thus we may 
restrict ourselves to an analysis of the cross section for the process in Fig. |^a, denoted as 
da a- To a very good approximation, this cross section can be expressed by the number of 
equivalent photons drti emitted from one of the nuclei, given by Eq. ([T^), and the cross 
section for the process 7^4 — > OM + A, denoted as cr^^, 

daa = drii a^A ■ (24) 

We thereby assume that the incident photon in the process OM + A is a virtual 

photon in the framework of the equivalent photon approximation and thus exhibits a small 
virtuality ^ < 4m^. Therefore, we neglect the virtuality of this incident photon in 
the cross section a^A- The subprocess •jA — > OM + A is described by the set of diagrams 
of Fig. ^. We calculate its cross section a^A in the region of large energies and relatively 
small transverse momenta |p^| of the produced OM. We have the kinematical conditions 
(the subscript 'th' denotes the threshold value) 

s-y = 2 gi • P2 ^ Sth = 4 m^M and |p^|~m^. 

We note, that a loop integral has to be evaluated for this subprocess. Its contribution is rather 
different from what would be obtained within the standard equivalent photon distribution 
for the two remaining photons. 

It is convenient to perform the calculations involved using the impact representation, 
which has been employed in QED and QCD for a number of processes with two photon 
or two gluon exchange (in the t-channel). More details on this approach are described in 
Refs. pl|-p^. In this representation the amplitude M^a which corresponds to the whole set 
of diagrams of Fig. ^ is written with an accuracy ~ m^/ in the form of a two-dimensional 
integral over the transverse components of the momentum of the virtual photon. 

The impact factors and Ja correspond to the upper and lower block of the diagrams 
in Fig. |[ The diagrams in Fig. |^ are regarded as being cut by the photon lines of the 
lower block, dividing the process into two partial virtual processes, 7 + 77 — OM (upper 
block) and A — 77 + A (lower block). The impact factor then corresponds to the virtual 
transition 7 + 77 — > OM, and Ja corresponds to the virtual transition in the lower block 
(A — » 77 + y4). The impact factor J a for a charged point-like particle was found in 0,1^ as 
Ja = 4:TraZ^. In our case we should take into account the shape of the nucleus and modify 
this impact factor according to 

Ja = AiraZ^ F{kl) F((p^ - k^)^) . (26) 
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The impact factor for the virtual transition 7 + 77 — > OM is similar to the impact 
factor for the virtual transition 7 + (7(7 — \1/ which was introduced for the description of the 
hard diffractive process 7g 1^^. Adjusting for the different couplings and masses, we 

immediately obtain 



m. 



m. 



ml + pf^ ml + (pix - kj 



(27) 



Here pi_L = p2± = 1/2 p±, and and eoM are the polarization vectors for the initial 
photon and the final state OM. iFiom Eq. (^) it follows that the helicity is conserved in 
the 7 —>■ OM transition. Therefore, the OM is transversely polarized and is produced in two 
polarization states only (not three states). 
We finally obtain the cross section as 



da. 



Z a 



<l>(pi)e^eoM 



2 (Pp± 



m. 



(28) 



where $(pi) is determined by an integral related to the amplitude M^a given in Eq. 
$(p5_) can be written as 



*(pi) = - 



(r + n)2p^ 



r2_i 



(r 



n 



;i + r)(l + rr2 



(29) 



where r is a two-dimensional vector with no physical dimension, over which the integration 
has to be performed, n is a unit vector defined by 



n 



and T is given by 



\P±\ 



4 ml 



After integrating over the azimuthal angle of OM, summing over the polarizations of the 
final state (OM spin states) and averaging over the polarizations of the initial state (photon 
polarizations), we obtain 



a, 



7A 



^TrZ^a^ A2 
B — 



m. 



m. 



(30) 



where the dimensionless constant B follows from 



B 



mA^yfdu 



(31) 



and A is the form factor scale defined in Eq. (^. 

The value of the constant B depends stronger on the shape of the form factor than the 
corresponding quantity for paradimuonium. We used a realistic form factor of the nuclear 
charge density p(r) for which we employed the model p5| 
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F(k2) = 7^ / d^r e^-"" p(r) 
Ze J 



with 



The parameters are 



Ze 1 , . 

^^"■^ = aT l+exp((r-i?)/a) ' 



= 1.18 fm , a = 0.53 fm . 



is the normahzation factor chosen such that / d^r p(r) = 1. 

The evaluation of the constant B is performed numerically on IBM RISC/6000 work- 
stations. Because the form factor Eq. (|32D is a function of with R oc A~^, the constant 
B has the same value for all nuclei considered in this paper. We obtain 

B = 0.85 . (33) 

It is useful to consider the sensitivity of the result on the choice of the form factor. With 
the approximate form factor given in Eq. (|^ the function $ is calculated in the Appendix 
analytically (with an additional approximation r = in the denominators of Eq. (pQ])). 
Further evaluation results in B = 0.93, this value is in fair agreement with the exact value 
from Eq. (^). 

Because the cross section of the subprocess •yA — > OM + A (cf. Eq. (|30D ) is energy 
independent in the discussed limit, the remaining integration of Eq. (^) is in fact an in- 
tegration over the equivalent photon spectrum only. Let uj2 be the total energy of both 
exchanged photons in Fig. ^. Then we have as in the previous section 4ujiuj2 = 4 (due to 

four-momentum conservation and the kinematics of the process) and 002 ~ A7 (due to the 
nuclear form factor). Thus a lower limit for Ui is 

2 2 

mf, mf, , , 

cji > ^ ^ . (34) 

The upper bound in this integration can be set to 00 due to the fast decrease of the equivalent 
photon spectrum at large energy. We obtain (using the notation 6 = m^/ (A7) introduced 
previously) 

Z'^a 

(Ta = H{S) a^A , (35) 

n 



where 



m2/(A7) 



7 d,, / \ roo rlT f 57 for RHIC, Au mode, 

H[5) = / ^ / ^ = / — f{x) = 202 for LHC, Pb mode, 
^-.L-. VA7; X 1 247 for LHC, Ca mode. 



(36) 



Finally, the cross section for the OM production is equal to (cf. Eqs. ( ^p3| , p0 



A' ... . . 

aoy, = 2^^ — BH[5). 37 
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The numerical values are 



0.021 for RHIC, Au mode , 
0.089 for LHC, Pb mode , 
0.000069 for LHC, Ca mode . 



(38) 



The ratio for the production cross section for the ortho and para states is given by 



Hence, we expect predominantly a production of para states in relativistic heavy ion colli- 
sions. 



The dimuonic atoms are neutral systems produced by a number of photons which are 
approximately on-shell and collinear to the colliding ions. So the angular spread of the 
dimuonia with respect to the beam-axis is of the order of 0{'y~^). The rapidity of the 
DM particles will be correspondingly high. Therefore the DM systems will not be observed 
directly by any detector with a low rapidity coverage. 

The dominant decay channels of DM are 77 (for PM) and e~^e~ (for OM, cf. Eq. (^). 
Dimuonium could be observed via detection of these decay products. We will investigate in 
this Section the influence of three sources of background on the prospective measurements, 

• photon and electron-positron background originating from hadronic processes in the 
interaction region, 

• free electron-positron pair production shadowing the signal from the decaying or- 
thodimuonium and 

• photon pair production by the two colliding nuclei shadowing the signal from the 
decaying paradimuonium. 

The above sources of background are expected to affect mainly the signal of those dimuonium 
atoms which decay in or near the interaction region of the heavy ion collision. 

First we consider photon and electron-positron background originating from inelastic 
hadronic processes in the interaction region. This background will affect the signal from 
both OM and PM atoms. In these inelastic processes, one or both of the nuclei dissolve to 
some extent. One can roughly divide these processes into two classes. The first are mainly 
hadronic processes, where the two nuclei collide and the strong interaction takes effect. The 
cross section for this class can be estimated as aAA ~ 4 A^^^app, which for the nuclei under 
investigation is in the range of five to seven barn. The second type is a photodissociation 
process caused by an energetic photon emitted from a larger distance by one of the nuclei. It 
induces nuclear reactions in the other nucleus on impact. The cross section of this photodis- 
sociation process depends crucially on the type of the colliding ions. One finds cross sections 
of roughly 85 barn for RHIC in the Au mode, 200 barn for LHC in the Pb mode and 3 




(39) 



IV. ESTIMATE OF THE BACKGROUND 



10 



barn for LHC in the Ca mode [|T^]. In the following we list the approximate luminosities per 



bunch crossing and the corresponding probability of hadronic events per bunch crossing 
'Ph = c"h -Cf), where the hadronic cross section ah is a combination of the purely hadronic and 
the photo-dissociation cross section. We obtain 



2.2 X 10^9 cm-2 , r 0.004 for RHIC, Au mode , 

Ch = i 3.75 X 10^0 cm-2 , Vt ^ I 0.075 for LHC, Pb mode , 



lO^-' cm"^ , [ 1 for LHC, Ca mode . 

By virtue of this figures we may conclude, that for Au and Pb the hadronic event rates are 
small enough to see DM-production in anticoincidences with the production of additional 
hadrons. In contrast this seems to be quite an impossible task for the Ca mode at LHC. 

A second, significant source of background for the decay of the OM is caused by the 
production of free e^e~ pairs by the two nuclei. In order to estimate this effect, we consider 
the e^e~ pair production via the standard two-photon mechanism. The cross section of this 
process, a^, is estimated with the well-known Racah-formula (see for details). We obtain 



(Te ~ 35 000 barn for RHIC and fXe ~ 225 000 barn for LHC. This is orders of magnitudes 
larger than the production cross section for the DM. A remedy for this problem might be 
a precise determination of the invariant mass of the electron-positron pair. The production 
cross section for an e+e~ pair having an invariant mass near 2m^ with mass spread Am 
is calculated using Eq. (|^) with the replacement of by the cross section for the process 
77 — >• e~^e~: 

AaiA, A, A, A, e- e') = G{6) fin - l") ^ ^ 1.8 x 10^^ a... (40) 

Because the orthodimuonium production rate is comparably low, even a realistic mass res- 
olution of 1 MeV would not fit our goal to distinguish the background from the signal of 
the orthodimuonia. Hence, it seems to be a very difficult task to observe the orthostate of 
dimuonium in heavy ion collisions. 

This situation is different for the parastate. The main non-hadronic background to the 
the deacy of the PM atoms is the two photon production process A1A2 —>■ A1A2 + jj. 
It can be described as the radiation of two (virtual) photons and subsequent light by light 
scattering (via an electronic loop). The cross section for this process is given by Eq. (|]) with 
the replacement of a^^ by the cross section of light by light scattering. It is five orders less 
than (Te and for invariant masses > 200 MeV of the photon pair it drops by five more orders. 
Nevertheless, the total free photon pair production cross section for all energies greater than 
200 MeV is still larger than the production cross section for PM by a factor of ~ a~^. To 
improve the situation one might again try to fix the invariant mass m as precisely as possible. 
In general, a relative precision Am/m at ^ 4m^ will lead to a signal-to-noise ratio of 
the order of Am/m. 

. / . . < . N Z^a^ Am 1 Am , . 

Aa A2 -> 7 7 = 0.95 G{b) ^ apM- 41 

2m^ m^ a m^ 

So the signal to background ratio becomes about 0.75 for mass resolution Am = 1 MeV. 
This corresponds to a determination of the invariant mass of the decay products of PM with 
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a precision of roughly 5 x 10~^. By contrast, for OM, a determination of the invariant mass 
of the electron-positron pair to an accuracy of the order of 10~^ would be necessary in order 
to reach a comparable signal-to-noise ratio. 

Another possibility to further improve the situation is to take into account only highly rel- 
ativistic dimuonium systems, which decay outside the nuclear collision region. The electron- 
positron or photon pair (for OM and PM, respectively) is produced outside the interaction 
region. DM atoms decay after travelling a typical decay length is £ ~ e/{2m^)cT (see Ta- 
ble H). This distance is increased for the excited dimuonia. The opportunity to observe 
excited dimuonium states in this approach will be better than that for the ground state de- 
spite the smaller production rate. The task left in this picture is to reconstruct the vertices 
of the decays, which necessitates a vertex detector. 



V. RESULTS AND DISCUSSION 



We have analyzed the production of bound states consisting of a muon and an antimuon 
in relativistic heavy ion collisions. The analysis of the parastate production was performed for 
the dominating two photon process. The two photon approximation describes the production 
of the parastate with an accuracy of (1 -r- 2)%. A novel three photon mechanism for the 
production of the orthostates was developed. The accuracy of this approximation is (6-i-12)%. 
The theoretical uncertainty of our results is primarily due to multiphoton processes. Other 
sources of uncertainty like the dependence on the nuclear form factor or corrections to the 
equivalent photon spectrum, have been analyzed in detail. They are on the level of one 
precent. Because multiphoton processes enhance the production rate, our results should be 
regarded lower bound on the total production. 

We obtained numerical results for the dimuonium production at the new heavy ion 
colliders RHIC and LHC The results for all colliders and for a set of atomic states are 
presented in Table |T|. In Table |l| we also consider the properties of the atomic 2P dimuonic 
state which is produced in atomic transitions from 3S and 4SQ 

The dimuonic atoms travel, after production, with small angular spread along the beam 
axis. Therefore they are detectable by their decay products only. In general the extraction 
of a signal from the experiment will be easier for paradimuonium than for orthodimuonium. 
The reasons are twofold, 

• the total production cross section is much larger and 

• the background is significantly reduced. 

As it has been shown in Sec. the photon pair background shadowing the paradimuo- 
nium signal is roughly five orders of magnitude smaller than the free electron-positron pair 
background shadowing the signal from orthodimuonium. Additionally, the total production 



The results are obtained according to the treatment of recoil effects in Ref. [26|. The atomic 
transitions in heavy fermionium has been discussed in Ref. [^^. The DM spectrum is considered 
in detail in Ref. [01. 
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cross sections for para states are larger than those for ortho states by a factor of 10 to 
100, depending on the coUider and the nucleus used (see Table |l|). We expect a favourable 
signal-to-noise ratio for the parastate if the energy of the photon pair can be determined 
with a precision of roughly 1 MeV. 
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APPENDIX A: PARADIMUONIUM LIFETIME 

We consider briefly the lifetime of the parastates of dimuonium. Because of the higher 
production rate of parastates, this is of interest in the context of possible experiments. The 
leading term is caused by the 77 decay, 

r(o)(,„is„) = i^. (Ai) 

Typical contributions to the next-to-leading order (NLO) corrections are depicted diagram- 
matically in Fig. ^. The corrections are evaluated in [|| as 

Ar^^°(ii^o) = 4.79 - r(°)(i^^o) (A2) 

TT 

and 

Ar^^°(2^5o) = 4.65 - r(°)(2^5o) . (A3) 

TT 

The next-to-next-to-leading order (NNLO) corrections include the large logarithmic factors 
ln(l/a;) and In^ (m^/me). We consider here these logarithmic terms. The ln(l/a) term is of 
the same form as for parapositronioum [p^ ], 

^pNNLo^^i^^) = 2aMn T^'^Kn'So) . (A4) 

The double mass ratio logarithm does not have an analogy in parapositronium. The relevant 
Feynman diagrams are depicted in Fig. 0. We obtain the result, 

AVr-in'So) = (1 + 2) In^ T^'Hn'So) . (A5) 

The terms (1+2) originate from the diagrams in Fig. ^ and Fig. 0b, respectively. The calcula- 
tion of the double logarithmic corrections is done here by the evaluation of the imaginary part 
of the diagrams in Fig. |^. The real part of these diagrams contributes to the hyperfine struc- 
ture in higher order. The final results for the lifetimes are r(l^S'o) = 0.59504(22) x lO^^^s 
and t{2^So) = 4.7619(17) x lO'^^ g_ 
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APPENDIX B: DEPENDENCE OF THE ORTHOSTATE PRODUCTION ON THE 

FORM FACTOR 



Employing the approximate form factor given in Eq. , we evaluate the integral Eq. (|3TD 
analytically. First we observe that the predominant contribution to B is caused by the region 
where p'j_/A.'^ < 1. In this region r = p\/{^Tn?^ < 0.1. Therefore we may put r = in 
Eq. (|9D. We integrate 

$(A^ = 1 / - ^) (Bi) 

^ ' 7r7 (r-n)2(r + n)2 [l + (M/4)(r-n)2][l + (M/4)(r + n)2] ■ ^ ' 



The relation 



'r + n)^ + (r - n)^ - 4 



proves to be useful for a simplification of the integrand. We can present the integral in 
Eq. ( pi|) in the form 

^{\^u) = -(2 + ai)/ii -2/22 + (4 + ai) Ji2. (B2) 

The integrals Uj are defined as 

If (fr 



TT J [(r + n)2 + ai] [(r - n)^ + aj 



where ai = A/u. The auxiliary parameter 02 = e ^ is introduced in order to regularize 
divergences in intermediate calculations. For the integrals lij we obtain after elementary 
integration 

1 , l + yiT^ ^ li 4 ^ 1 (4 + ai)2 
ill = , In = , /22 = - In - , J12 = — In 



vT+oi ^/al 2 e 4 + ai aie 

Substituting these expressions into Eq. ( p2D we obtain the result 



u) = In il±^ - 4±2^ In , (B3) 

in which the dependence on the auxiliary parameter e vanishes. Using this method we obtain 
a result of 5 = 0.93 in contrast to B = 0.85 with the exact form factor (cf. Eq. (|33D ). 
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TABLES 



Ooliider 


Nucleus 


Nuclear 


Luminosity 


Lorentz 


Bunch 






charge 


C 


factor 


length 






{Z) 


[cm~^s~-^] 


(7) 


[cm] 


RHIC 


AuAu 


79 


2 X 10^6 


108 


12 


LHC 


PbPb 


82 


3 X 1027 


2980 


7.5 


LHC 


CaCa 


20 


4 X 10^° 


3750 


7.5 



TABLE L Experimental parameters of RHIC and LHC which must be taken into account 
for the production of dimuonium. The bunch length of 7.5 cm for the CaCa-channel at LHC is an 
estimate. 



ATOMIC STATE PROPERTIES 



ESTIMATED PRODUCTION PER YEAR 



Atom 


State 


jPC 


cr 


Decay 


RHIC 


LHC 


LHC 








[cm] 


mode 


Au-Au 


Pb-Pb 


Ca-Ca 


PM 


l^^o 


0-+ 


0.0178 


77 


310 


40000 


260000 


PM 


2^5o 


0-+ 


0.143 


77 


40 


5000 


33000 


PM 


3^5o 


0-+ 


0.483 


77 


12 


1500 


9800 


PM 


4^So 


0-+ 


1.14 


77 


5 


630 


4100 


PM 


8^So 


0-+ 


9 


77 




79 


520 


PM 


lO^^o 


0-+ 


18 


77 




40 


260 


OM 


13^1 


1— 


0.0538 




43 


2700 


2800 


OM 


2^51 


1— 


0.430 




5 


330 


340 


OM 


3^51 


1— 


1.45 


e^e~ 


2 


100 


100 


PM 


2^ Pi 


1+- 


0.462 


l'5o7 




60 


400 



TABLE II. Main properties of atomic states of dimuonium and their estimated production at 
LHC and RHIC per year (running time per year in our calculation is 10^ s). The decay mode given 
here is the dominant mode which is most important for the detection. 
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FIGURES 




FIG. 1. Diagrams for two and three photon production mechanisms of fermion pairs in rel- 
ativistic heavy ion coUisions. In the case of dimuonium, the fermion pair is produced in a bound 
state. 



=1 — 




I 



FIG. 2. Two photon production of paradimuonium by relativistic heavy nuclei. 



FIG. 3. Distribution of paradimuonium produced at LHC in the Pb mode over the energy (in 
GeV). The distribution is normaUzed to the annual production rate of 40000 particles. The median 
of the distribution is at 1.12 GeV. 
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diagrams originating 
I from the symmetriza- 
~ tion of the photon 
wave functions. 



fqi-k 




+ 



diagrams originating 
from the symmetriza- 
tion of the photon 
wave functions. 



(b) 



FIG. 4. Orthodimuonium production by a three photon fusion process. 




FIG. 5. Orthodimuonium production by a three photon fusion process. P2 denotes the nuclear 
momentum, p is the momentum of the dimuonium system. 
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